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[57] ABSTRACT 

A method is provided for determining surface shear 
magnitude and direction at every point on a surface. 
The surface is covered with a shear stress-sensitive 
liquid crystal coating and illuminated by white light 
from a normal direction. A video camera is positioned 
at an oblique angle above the surface to observe the 
color of the liquid crystal at that angle. The shear mag- 
nitude and direction are derived from the color infor- 
mation. A method of calibrating the device is also pro- 
vided. 

9 Claims, 5 Drawing Sheets 
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METHOD FOR MEASURING SURFACE SHEAR 
STRESS MAGNITUDE AND DIRECTION USING 
LIQUID CRYSTAL COATINGS 

5 

ORIGIN OF THE INVENTION 

The invention described herein was made by an em- 
ployee of the United States Government and may be 
manufactured and used by or for the Government for 
governmental purposes without the payment of any 1 
royalties thereon or therefor. 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This is a Continuation-in-Part of application Ser. No. 15 
08/031,972, filed Mar. 16, 1993, now U.S. Pat. No. 
5,394,752. 

FIELD OF THE INVENTION 

20 

The present invention is directed to aerodynamic 
testing using liquid crystals, and, more particularly, is a 
method for obtaining full-field measurements of both 
the magnitude and direction of surface shear stress in 
such testing. 25 

BACKGROUND OF THE INVENTION 

In fluid mechanics and aerodynamics research, much 
valuable information can be gained from visualizations 
of both the outer flow and surface shear stress patterns 3Q 
over solid bodies immersed in fluid streams. The liquid 
crystal coating method, i.e., a method wherein a liquid 
crystal coating is applied to a surface under test, is a 
diagnostic technique that can provide areal visualiza- 
tions of instantaneous shear stress distributions on sur- 35 
faces in dynamic flow fields with a response that is 
rapid, continuous and reversible. Reference is made, for 
example, to Klein, E. J., “Liquid Crystals in Aerody- 
namic Testing,” Astronautics and Aeronautics, Vol. 6, 
July 1968, pp. 70-73. 40 

Liquid crystals are highly anisotropic “fluids” that 
exist between the solid and isotropic liquid phases of 
some organic compounds. As such, these crystals ex- 
hibit optical properties characteristic of a crystalline 
(solid) state, while displaying mechanical properties 45 
characteristic of a liquid state. Typically, in flow- visual- 
ization applications, a mixture of one part liquid crystals 
to nine parts solvent (presently, Freon) is sprayed on the 
aerodynamic surface under study. A smooth, flat-black 
surface is essential for color contrast and must be kept 50 
free of grease and other chemical contaminants. Rec- 
ommended applications (after spray losses) are about 10 
to 20 ml liquid crystals, measured prior to mixing with 
the solvent, to each square meter of surface area. The 
solvent evaporates, leaving a uniform thin film of liquid 55 
crystals whose thickness, based on mass conservation 
and estimated spray losses, is approximately 10 to 20ju, 
(0.0004 to 0.0008 inch). Once aligned by shear, the mol- 
ecules within the liquid crystal coating selectively scat- 
ter incident white light as a spectrum of discrete colors, 60 
with each color at a discrete angle (orientation) relative 
to the surface. For “thermochromic” liquid crystal 
compounds, this molecular structure, and thus the light 
scattering capability of the coating, responds to both 
temperature and shear stress. For newly-formulated, 65 
shear-stress-sensitive/temperature-insensitive com- 
pounds, such as Hallcrest BCN/192 and CN/R3, “color 
play ” (i.e., discerned color changes at a fixed angle of 
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observation, for a fixed angle of illumination) results 
solely from the application of shear stress. 

While it is now known that this technique can be 
calibrated (under carefully controlled conditions) to 
measure surface shear stress magnitudes, two important 
issues remain: time response and “directional sensitivi- 
ties,” i.e., sensitivity to illumination and viewing angles, 
as well as to the instantaneous shear stress direction. 

The time-response issue has been investigated (see 
Parmar, D. S., “A Novel Technique for Response 
Function Determination of Shear Sensitive Cholesteric 
Liquid Crystals for Boundary Layer Investigations,” 
Review of Scientific Instruments, Vol. 62, No. 6, June 
1991, pp. 1596-1608) by placing a liquid crystal layer 
(about 100 p, thick) between two optical glass plates and 
applying known and transient shear forces via a dis- 
placement of one plate. Liquid crystal delay, rise and 
relaxation time constants were measured as a function 
of the monochromatic wavelength of the incident light. 
Time constants in the range of 10 to 100 milliseconds 
were generally observed, with minimum values being 
on the order of 3 milliseconds. The extrapolation (or 
applicability) of these results to actual fluid mechanic 
applications, wherein liquid crystal coating thicknesses 
are an order of magnitude less, remains an open ques- 
tion. 

Another approach to characterizing the time- 
response and/or flow-direction-indication capabilities 
of the liquid crystal technique is to expose the coating to 
transient viscous flows of known time scales and/or 
known shear directions. Oscillating airfoil experiments 
provided some initial results concerning both issues (see 
Reda, D. C., “Liquid Crystals for Unsteady Surface 
Shear Stress Visualization,” AIAA Paper 88-3841, July 
1988 and Reda, D. C., “Observations of Dynamic Stall 
Phenomena Using Liquid Crystal Coatings,” AIAA 
Journal, Vol 29, No. 2, February 1991, pp. 308-310). 

In general, while previous investigators have utilized 
the liquid crystal technique, the focus of such research 
has been on attempting to determine parameters related 
to the magnitude of shear stress, for example the sudden 
increase in shear stress levels at the transition from 
laminar to turbulent flow. Reference is made, for exam- 
ple, to U.S. Pat. No. 4,774,835 (Holmes) and to Holmes 
et al, “A New Method for Laminar Boundary Layer 
Transition Visualization in Flight: Color Changes in 
Liquid Crystal Coatings,” NASA Technical Memoran- 
dum, 7666, January 1986. 

Currently, numerous point measurement techniques 
are available to determine surface shear stress magni- 
tude. These techniques only allow discrete point mea- 
surements of surface shear stress magnitude. Several of 
these methods can also be used to measure local shear 
stress direction. However, no full-surface vector mea- 
surement capability presently exists. 

SUMMARY OF THE INVENTION 

It is an object of the present invention to use the 
liquid crystal coating (LCC) method as a quantitative 
technique for measuring surface shear stress vector 
fields in three-dimensional flows. 

The present invention is based on the discovery or 
inventive appreciation that when the liquid crystal coat- 
ing is subjected to shear (e.g., by a wind tunnel flow 
over the surface of the liquid crystals), the resulting 
color-play response is dependent on both the magnitude 
and the direction of the shear. By directing a beam of 
white light onto an LCC surface to cause the white light 
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to be dispersed from the surface, generating a calibra- 
tion curve of color versus absolute shear magnitude, 
rotating an observer from a setting used for the calibra- 
tion step, and obtaining a plurality of images by the 
observer along a circumferential arc, the shear direction 5 
may be determined from the plurality of images and, 
from the color corresponding to the shear direction, the 
shear magnitude may be determined from the calibra- 
tion curve. 

Advantageously, in generating the calibration curve, 10 
the maximum shear level is applied to the surface, the 
observer is rotated about the surface until a maximum 
color change signal is located, an above plane view 
angle is adjusted until a further maximum color change 
signal is located. The calibration step then includes 15 
generating a calibration curve for a plurality of para- 
metric values of shear magnitude. 

Advantageously, the shear direction may be deter- 
mined by plotting color versus circumferential view 
angle, curve fitting a symmetric function along the plot, 20 
and determining, from the minimum of said symmetric 
function, the shear direction and corresponding color. 

Other features and advantages of the invention will 
be se forth in, or apparent from, the following detailed 
description. 25 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic representation of the dispersed 
(reflected) light spectrum produced by directing a beam 
of white light obliquely onto a surface coated with an 30 
aligned liquid crystal layer; 

FIG. 2(a) and 2(b) are views corresponding to that of 
FIG. 1 showing the rotation of the reflected light spec- 
trum in response to shear stresses imparted to the liquid 
crystal coating in respective opposite directions; 35 

FIG. 3 is a top plan view of a measuring arrangement 
of the present invention; 

FIG. 4 is a side elevational view, partially in cross 
section, of the measuring arrangement shown in FIG. 3; 

FIG. 5 is a CIE chromaticity diagram for numerous 40 
shear magnitudes at /3=0°; 

FIG. 6 is a plot of dominant wavelength versus rela- 
tive surface shear stress magnitude for limiting circum- 
ferential view angles of fi — 0° and 180°; 

FIG. 7 is a plot of dominant wavelength versus rela- 45 
tive circumferential view angle between observer and 
shear vector, with relative surface shear stress magni- 
tude as the parameter; and 

FIGS. outline the four steps necessary to mea- 

sure surface shear stress vector fields using the liquid 50 
crystal coating technique. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

Referring to FIG. 1, a schematic representation is 55 
provided of the “color play” provided by a liquid crys- 
tal coating. In FIG. 1, an aligned liquid crystal coating 
10 is provided on a black surface 12 beneath coating 10, 
and a source 14 produces a beam of white light indi- 
cated at 16. The reflected spectrum 18 from coating 10 60 
comprises, as discussed above, a series of discrete colors 
or color bands. More specifically, in the schematic rep- 
resentation provided in FIG. 1, these bands comprise a 
blue (turquoise) band 20, a green band 22, a yellow band 
24 and a red band 26, and these bands are produced in 65 
addition to a reflected glare line indicated 28. An ob- 
server, in the form of a recording device 30 comprising 
a movie and/or video camera, is positioned near “mid- 
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spectrum,” i.e. just before the yellow-to-green transi- 
tion angle so that the observer sees the yellow color. In 
this simplest of all possible arrangements, the light, 
camera and flow direction are all in the same plane. The 
technical paper “Time and Flow-Direction Responses 
of Shear-Stress-Sensitive Liquid Crystal Coatings,” D. 
C. Reda, J. J. Muratore, Jr., and J. T. Heineck, AIAA 
Journal, v. 32, April 1994, pp. 693-700, which is hereby 
incorporated by reference, describes calibration proce- 
dures for more complex experimental arrangements. 

Referring to FIGS. 2(a) and 2(6), there are shown 
schematic representations similar to that of FIG. 1, 
illustrating the spectrum rotations produced by the 
application of shear stresses. In particular, FIG. 2(a) 
depicts a situation wherein shear, i.e., shear stress, indi- 
cated by arrow SI, is exerted from left to right and as 
illustrated, under these circumstances the reflected 
spectrum 18 rotates counterclockwise so that the blue 
band 20 is viewed or observed by camera 30 instead of 
the yellow band as in FIG. 1. In contrast, 2(b) illustrates 
a situation wherein shear, indicated by arrow S2, is 
exerted from right to left and illustrates the fact that the 
spectrum 18 rotates clockwise and the red band 26 is 
viewed by recording device 30. Thus, it will be appreci- 
ated from the foregoing that by determining the color 
band observed by the observer (recording device) 30, a 
determination can be made of the direction of the shear 
applied to the liquid crystal coating 10. 

As demonstrated above, liquid crystal coating (LCC) 
color change response to shear depends on the shear 
stress direction. It is further known that LCC color play 
response depends on the shear stress magnitude. Based 
on this knowledge that the color-play response depends 
on both the shear stress magnitude and direction, the 
arrangement shown in FIGS. 3 and 4 allows the LCC 
method to be converted from a flow visualization tool 
to a quantitative measurement technique for surface 
shear stress vectors. The technical paper “Measurement 
of Surface Shear Stress Vectors Using Liquid Crystal 
Coatings,” D. C. Reda and J. J. Muratore, Jr. AIAA , 
Journal, v. 32, August 1994, pp. 1576-1582, which is 
hereby incorporated by reference, describes the mea- 
surement methodology. 

In contrast to the arrangement shown in FIGS. 
1-2(6) in which oblique illumination is used, FIGS. 3 
and 4 illustrate the use of normal illumination. Under 
normal illumination, any point exposed to a shear vec- 
tor with a component directed away from the observer 
exhibits a color play response, i.e. a shift from the no 
shear color (orange) towards the blue end of the visible 
spectrum. Conversely, any point exposed to a shear 
vector with a component directed toward the observer 
exhibits a non-color play response (always character- 
ized by a rusty brown appearance). Further, this two- 
lobe color response rotates with the observer, i.e., it is 
identical for all circumferential view angles. 

For the arrangement in FIGS. 3 and 4, a liquid crystal 
mixture having a viscosity and full-spectrum color-play 
response consistent with the absolute shear levels to be 
encountered must be selected and a coating 78 applied, 
here to a plate 80. Plate 80 is supported by a support 
stand 82, which in turn is mounted on the optics table 
86, The surface of plate 80, assumed planar for the sake 
of this explanation, is then illuminated with white light 
62 from white light source 60 in the normal direction 
and the maximum shear level is applied to the surface. 

A pair of pressure regulators 70, 72, connected in 
series, control the pressure difference used to drive the 
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flow through the circular jet 68. The value of this pres- 
sure difference is indicated by a calibrated gauge 76. A 
valve 74 is used to supply the pressure difference. 

In order to calibrate the liquid crystal color play 
response to shear, surface shear vectors of known direc- 
tion (relative to the observer) and known relative mag- 
nitude are imposed upon coated surface 80 via jet 68. In 
such flows, the skin friction coefficient, Cf^T/O.SpVj 2 , 
is essentially constant in the near exit region over a wide 
range of Reynolds numbers. Thus, for a constant den- 
sity (p), the magnitude of the surface shear stress vector 
(r) on the jet centerline scales linearly with the square 
of the jet exit velocity (V/ 2 ) which in turn scales linearly 
with the pressure difference used to drive the flow. 
Relative changes in jet pressure, measured by gauge 76, 
define the relative changes in shear magnitude (r/r r ), 
where the subscript r denotes some arbitrary “lowest 
shear” reference level. 

An observer, here a fiber optic probe 64, is positioned 
above the test surface at a nominal above-plane 
(oblique) view angle (indicated at a p ), e.g., 0^=30°, 
facing downstream relative to the nominal freestream 
velocity vector. The measurement point of the fiber 
optic probe 64 is on the centerline of the wall jet flow 
and this point is always coincident with the center of 
rotation of the probe-traversing platform, thereby al- 
lowing the relative view angle between the shear vector 
and the observer ( fi ) to be systematically varied. Light 
captured by probe 64 is input to a spectrophotometer 66 
which disperses the received light into its spectral con- 
tent. Measurements of irradiance from the coating 78, as 
a function of the wavelength of light, are thus obtained 
over the entire visible spectrum. 

As shear magnitude is increased, the peak intensity of 
the scattered light decreases and, more importantly, the 
wavelength corresponding to this peak intensity shifts 
to lower wavelengths (indicating a shift in color from 
orange, through yellow and green to blue). Using estab- 
lished methods, such irradiance measurements are con- 
verted into standard chromaticity coordinates (x,y) 
which quantify the color of the coating at each test 
condition. Graphical determination of dominant wave- 
lengths Xd is shown in FIG. 5, which is a CIE chroma- 
ticity diagram. FIG. 5 clearly shows the relationship 
between the color of the liquid crystal coating and the 
relative magnitude of shear stress applied at /3=0°. 
Dominant wavelengths (Xd), obtained by extrapolation 
from white light (W) through the measured chromatic- 
ity coordinates (the solid dots in FIG. 5) to the extremes 
of the chromaticity diagram, have been found to be 
entirely consistent with the observed peak-power wave- 
lengths and the visually observed colors. A single cali- 
bration curve of color (dominant wavelength \£>or hue) 
versus relative shear stress magnitude can then be ob- 
tained, as shown in FIG. 6; conventional point-measure- 
ment techniques should be used to measure absolute 
shear magnitude r in actual experiments. 

FIG. 6 shows a plot of measured minimum dominant 
wavelength ( XDmin ) versus relative shear stress magni- 
tude for the situation where the observer and the vector 
are perfectly aligned and flow is away from the ob- 
server (at /3=0°). Under this condition, color change is 
seen to scale linearly with increasing shear stress magni- 
tude for l^r/r r ^8. At higher shear rates, the color- 
play response becomes nonlinear. 

For non-planar test surfaces, a family of such color 
versus shear magnitude calibration curves would have 
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to be generated to parametrically cover those effective 
(above-tangent-plane) view angles encountered. 

FIG. 7 shows plots of dominant wavelength (Xd) 
versus the relative, circumferential-view angle between 
5 the shear vector and the observer (/3, where positive 
angles are measured clockwise from the local line of 
sight). The probe angle, relative to the plane of the test 
surface (a^) is held constant. The parameter in this fam- 
ily of curves is the relative magnitude of the shear vec- 
10 tor (r/r r ). The key for the plot is as follows: O 
(r/T r )= 1; ■ (T/r r )-2; A (r/r r ) = 4; #( t / t >)=6; and □ 
(r/r r )= 8. Two principal results are evident: (1) for a 
given shear stress magnitude, the minimum dominant 
wavelength (i.e., the maximum color change) is mea- 
15 sured when the shear vector is aligned with, and di- 
rected away from, the observer (at /3=0°); (2) for a 
given shear stress magnitude, changes in relative cir- 
cumferential-view angle /3 to either side of the vector- 
/observer aligned position result in symmetric increases 
20 in measured dominant wavelength. While the solid lines 
in FIG. 7 indicate Gaussian curve fits, which fit the 
measured data sets well for —90° ^(3^ +90°, it is noted 
that the data may also be fit by a simple second-order 
polynomial curve over the limited range 
25 — 60°=i/3= +60°. 

Influences of point-to-point variations for a single 
coating as well as coating-to-coating variations are not 
great. At relative shear magnitudes of 4, 6, and 8, mea- 
surements of dominant wavelength were repeated to 
30 within order 2 to 3 nm over the entire (3 range shown in 
FIG. 7. Similar small variations in Xd, of order 2 nm or 
less, were measured after six-minute exposures to the 
highest shear levels of (t/t>)=8. 

A final element to be addressed concerns the ap- 
35 proach by which the method of the present invention 
can be expanded from a point-measurement technique 
to a full-surface, or areal, technique. To achieve this 
capability, the fiber optic probe and spectrophotometer 
are replaced by a three-chip, co-site sampled, RGB 
40 digital color camera (denoted by subscript c) linked to a 
full-image frame grabber and a supporting computer. 
Calibration and data acquisition are accomplished in 
exactly the same manner as used in the point-measure- 
ment techniques. FIGS. Sa-Sd schematically outline the 
45 four steps required to measure shear vector orientation 
(4> r ) and magnitude (r) at every point on a test surface. 

In the first step, the calibration curve for (3=0°, as 
shown in FIG. 8 u, is generated for a given a c , for the 
calibration vector aligned with and directed away from 
50 the camera (<f> c =4> r ). In the second step, three or more 
full surface images are acquired, maintaining a constant 
a Cf at known <J> C angles along a circumferential arc en- 
compassing all possible vector directions encountered 
in the three-dimensional, steady state flow under inves- 
55 tigation, as shown in FIG. 8 A 

The three RGB voltages for every pixel in each full- 
surface image are converted into a color value consis- 
tent with the conversion used in FIG. 8(c) In the third 
step, color versus circumferential angle <f> is plotted and 
60 a curve fit is obtained as shown in FIG. 8c, for each 
point on the surface. The correspondence between im- 
age-pixel locations and surface-point locations must be 
established for all view-angle combinations used to 
record the images. The curve fit must be a function that 
65 is symmetric about the minimum color value, e.g., 
Gaussian or parabolic. The angle at which the minimum 
color (or maximum color change) occurs defines the 
vector direction <f > T . 



5,438,879 


7 

In the fourth step, for each physical point on the 
surface, the curve-fit-defined minimum color value de- 
termined in the third step is input to the color versus 
shear magnitude calibration curve of the first step, as 
shown in FIG. 8 d, in order to define the vector magni- 5 
tude r. 

The test methodology outlined above does not re- 
quire calibration curves at any 0 orientation other than 
the observer/ vector aligned position (at 0=0°); it only 
requires the knowledge that the minimum dominant \q 
wavelength for any given shear stress magnitude occurs 
when the shear vector is aligned with, and directed 
away from, the observer. In other words, the family of 
curves represented by FIG. 7 is not required to apply 
the present measurement technique. Extension of this 15 
technique to unsteady (time-varying) flows would re- 
quire fixed/synchronized cameras to be positioned at 
three or more circumferential-view angles encompass- 
ing the vector directions to be measured. 

In summary, an LCC surface is illuminated with 20 
white light from the normal direction and a maximum 
shear level is applied to the surface. An observer is 
rotated about the surface until a maximum color change 
is located. At this orientation, the above plane view 
angle is adjusted in order to maximize the color change 2 „ 
signal. A color versus absolute shear stress magnitude 
calibration curve is obtained. Then, while the above 
plane view angle is maintained, images at known cir- 
cumferential view angles are obtained. These images 
are used to plot color versus circumferential view angle 
at each surface point and a curve fit is used to locate the 3 
minimum color value. The angle at which this minimum 
color value (or maximum color change) occurs is the 
vector orientation at each surface point. The calibration 
curve is then used to determine the magnitude of the 
shear vector corresponding to the minimum color value 35 
at each surface point. 

Resolution of the measurement technique, to ±1° in 
orientations and ±5% in magnitudes, is best achieved 
using five or more full-surface images and Gaussian 
curve fits. 40 

Alternatively, it is mathematically possible to deter- 
mine the two unknowns (magnitude and direction) from 
just two measurements of the coating color-play re- 
sponse to shear. Such a method, however, requires 
considerably more effort in the calibration process. 45 
Rather than generating a single calibration curve of 
color versus shear magnitude at 0— 0°, the entire family 
of color versus 0 curves for many parametric values of 
shear magnitude must be obtained, as shown in the plot 
of FIG. 7. 50 

Two fixed/synchronized cameras (not shown) would 
be deployed at a c , <f> c i, and a c , <j> C 2 such that the two 
lines-of-sight encompassed the vector directions to be 
measured. The in-plane angular spacing A<j> c between 
the two cameras would thus be known. 55 

Converting the two color measurements (Xjd\ and 
Xm) into a vector magnitude would then involve inter- 
polation within the family of curves depicted by FIG. 7. 
The measured vector magnitude would be defined 
when the two measured colors were found to lie on the 60 
same shear-magnitude curve subject to the condition 
|/*i | + 1 02 1 =A<j> C) where 1 02 [ < 1 01 1 when X D2 <Xm 
or |0i | <|02 1 when \m<\m , 

By experimental arrangement, one value of 0 would 
be positive and one would be negative. Having interpo- 65 
lated within the family of calibration curves to define 
the correct shear magnitude, the two values of 0 could 
then be read directly from the two X/> intercepts with 
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this “correct-magnitude” curve. The measured vector 
orientation would then be determined from geometry 
using either pair of angles, (j>ci and 0 i or <j>c2 and 02- 

It will be understood by those skilled in the art that 
although the invention has been described relative to 
exemplary preferred embodiments thereof, variations 
and modifications can be effected in these embodiments 
without departing from the scope and spirit of the in- 
vention. 

What is claimed is: 

1. A method for determining surface shear stress 
direction and magnitude at every point on a surface, 
said method comprising the following steps: 

providing a liquid crystal coating on a surface; 

directing a beam of white light perpendicular to said 
surface onto said surface; 

positioning a color-measuring camera at an oblique- 
view angle above said surface; 

generating, in a calibration step, a calibration curve of 
color versus absolute shear magnitude for calibra- 
tion shear vectors aligned with and directed away 
from said camera; 

rotating said camera from the setting used for said 
calibration step; 

obtaining a plurality of images by said camera along a 
circumferential arc; 

determining said shear direction from said plurality of 
images; and 

determining, from the color corresponding to said 
shear directions, said shear magnitudes from said 
calibration curve. 

2. A method as claimed in claim 1, wherein said cali- 
bration step comprises applying a maximum shear level 
to said surface. 

3. A method as claimed in claim 1, wherein said cali- 
bration step comprises rotating said camera about said 
surface, said oblique view angle remaining fixed, until a 
maximum color change signal is located. 

4. A method as claimed in claim 3, wherein said cali- 
bration step further comprises adjusting the above- 
plane view angle until a further maximum color change 
signal is located, the resulting above-plane view angle 
to be maintained constant for both calibration and data 
acquisition. 

5. A method as claimed in claim 1 wherein said cali- 
bration step comprises using point-measurement tech- 
niques to measure absolute shear stress magnitudes. 

6. A method as claimed in claim 1, wherein said plu- 
rality of images for data acquisition includes three or 
more images. 

7. A method as claimed in claim 1, wherein the step of 
determining shear direction at each physical point on 
the surface comprises plotting color versus circumfer- 
ential view angle, curve fitting a symmetric function 
along said plot, and determining, from the minimum of 
said symmetric function, said shear direction and corre- 
sponding color. 

8. A method as claimed in claim 7, wherein said plu- 
rality of images are full-surface images and further com- 
prising plotting color versus circumferential view angle 
for identical surface points taken from each of said full- 
surface images. 

9. A method as claimed in claim 7, wherein said cor- 
responding color at said minimum of said symmetric 
function can be used in conjunction with said calibra- 
tion curve to determine the shear magnitude at each 
physical point recorded in each of said full-surface im- 
ages. 


* Sjc * * * 



